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First-principles electronic structure calculations have been performed for the double perovskite 
Bi2CoMn06 in its non-centrosymmetric polar state using generalized gradient approximation plus 
Hubbard U approach. We find that while Co is in a high spin state, Mn is in a intermediate spin 
state. The calculated dynamical charge tensors are anisotropic reflecting a low-symmetry structure 
of the compound. Magnetic structure dependent phonon frequencies indicate the presence of spin- 
phonon coupling. Using Berry phase method, we obtain a spontaneous electronic polarization of 
5.88 fiC/cm^ which is close to the experimental value observed for a similar compound, Bi2NiMn06. 

PACS numbers: 75.50.Pp,75.70.-i,71.70.Gm 



Magnetic double perovskites with formula A2BB'06 
(A is an alkaline rare-earth, B and B' are transition met- 
als) have attracted a lot of attention in the last decade 
due to their potential in spintronic applications specially 
since the discovery of low field magneto-resistance in 
Sr2FeMo06 and Sr2FeRe06.'^. Depending on the sto- 
ichiometry and temperature, a number of double per- 
ovskites show inter-related magnetic, structural and elec- 
tronic phases. Although, many of these materials or- 
der ferromagnetically at high temperature, very few of 
them show ferroelectricitySHZl, A perpetual interest exists 
to identify materials having multiferroic properties with 
the simultaneous existence of ferroelectric and magnetic 
ordeil^. Here we have studied Bi2CoMn06 that crystal- 
lizes in the double perovskite structure. It is among a 
few, rare materials which show existence of both ferro- 
electric (FE) as well as ferromagnetic (FM) behavior. 
Due to the presence of 6s^ lone pair of Bi'^^, there is 
a non-centrosymmetric (space group C2 (no. 5)) struc- 
tural distortion which leads to a FE transition temper- 
ature around 600 kP, above which the centrosymmetric 
structure appears with space group C2/TO(no. 12). Be- 
low the magnetic ordering temperature of 95 K, ferro- 
magnetism appears due to the 180° cation-anion-cation 
ferromagnetic superexchange interaction^*^ in the chain 
Co"*"-0^-Mn+ where Co and Mn ions are ordered in 
a rock-salt configuration. A recent experimental study 
has claimed a strong dependence of Tc on the mon- 
oclinic angles] ^j^g present study, we concentrate 
on the C2 structure of Bi2CoMn06 (BCMO) as first- 
principles calculations have not been reported in the liter- 
ature for either of the two structural phases. Calculations 
were performed using density functional theory (DFT) 
within the Perdew-Burke-Ernzerhof generalized gradient 
approximation (GGA) on the basis of the projector aug- 
mented wave method implemented in the Vienna ab ini- 
tio simulation package (VASPJP^. The wave function was 
expanded in plane waves using 500 eV energy cut-off. For 
Bi and Mn, we have used pseudopotentials with semi-core 
electrons. Brillouin zone integrations were carried out us- 
ing a 4 X 6 X 4 Monkhorst-Pack fc-point mesh. To treat 
strong electron interaction effect in Co and Mn we have 




FIG. 1: (Color online) Crystal structure of BCMO, as viewed 
along the short-axis b, showing octahedral tilt patterns of 
CoOe (blue) and MnOe (red) octahedra along with Bi as pur- 
ple balls. 



carried out GGA-KU calculation in the Hubbard model 
framework with the choice of Coulomb parameter U=3 
eV for Co-d and U=4 eV for Mn-d orbitals. The value 
of the exchange parameter J used was 1 eV for both the 
transition metals. The unit cell consists of four formula 
units in the C2 phase. In Wyckoff notations, the four Co 
atoms occupy 2a and 2b sites respectively while Bi, Mn 
and O occupy 4c sites. In Glazer's notation, the tilt pat- 
tern of the octahedra is given by a~b"c~ (Fig. 1) which 
imply that there is no tilting of octahedra along [010] 
direction, while the neighboring octahedra tilt in oppo- 
site directions to each other along both [100] and [001] 
directions. We have performed the structural optimiza- 
tion involving the relaxation of ionic positions as well as 
cell shape and volume. Forces were relaxed up to values 
smaller than 0.005 eV/A. For the calculation of density 
of states we used a fc-point mesh of 6 x 6 x 6 with 54 
k-points in the irreducible Brillouin zone. The dielectric 
and dynamical properties such as Born effective charge 
(BEC) matrices and zon-center phonon frequencies were 
calculated using density functional perturbation theory 
(DFPT)i21. 
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It should be noted that we have used GGA+U method 
in the present study as GGA calculations are unable to 
yield an insulating state. Our optimized lattice param- 
eters are a= 9.62 A, b= 5.50 A and c= 9.74 A(unit cell 
volume of 490.21 A^). The corresponding experimental 
values are 9.49 A, 5.45 A and 9.57 A respectivelj^^l (unit 
cell volume of 471.13 A'^). The calculated monoclinic 
angle /3 between vectors a and c, is 108.15°, which is 
slightly overestimated with respect to the experimental 
one (107.60°). 
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FIG. 2: (Color online) Mn and Co DOSs projected on d- 
orbitals for (a) FM and (b) FiM configurations. The hori- 
zontal blue lines indicate separation between the occupied Co 
and unoccupied Mn states in the minority spin. See text for 
details. 

In Fig. 2, we show the density of states (DOS) obtained 
by GGA-I-U method for the ferromagnetic and ferrimag- 
netic (FiM) states. The choice of our U values opens up 
a band gap of around 0.3 cV in the FM state. The top 
of the valence band is dominated by Co-d and 0-p states 
(not shown) while the Mn states are pushed about 1 eV 
below the Fermi energy. On the other hand the bottom of 
the conduction band is due to Mn-d and 0-2p states. A 
comparison of the total energies of iso-structural FM and 
FiM structures yields that the FM state is 163 meV/f.u. 
lower in energy compared to the FiM state. The mag- 
netic moment of Co (both in 2a and 2b sites) is found to 
be 2.5 \iB^ while Mn carries a magnetic moment of 3.2 
PkBi yielding a magnetic moment of 6 [IbI'^-^- in the FM 
state. The analysis of orbital projected charges in both 
FM and FiM magnetic structures shows that Co is in a 
high spin state with around 7 d electrons (Co^+ state), 
while the Mn sphere contains around 5 d electrons. It 
seems that Mn is in an oxidation state which is associ- 
ated with an ionicity of about 2+ rather than 4+, which 
is the expected ionic state from the counting of oxidation 
states. This deviation from the expected ionicity may 



TABLE L Born effective charges of Bi, Co, Mn and O (one 
in the Co-O-Mn chain) for FM (left) and FiM (right) config- 
urations. 



be attributed to the low symmetry structure and small 
Mn-0 bond lengths. Similar ill-defined oxidation states 
of Mn and O were mentioned by Ciucivara el. aP for 
an iso-structural compound Bi2NiMn06. The structural 
analysis shows that in the FM state, LMn — O — Co(2a) 
is equal to 159.45° while the angle LMn - O - Co(26) 
is 150.83° giving an average LMn — O — Co angle of 
155.14°. In both the Mn-O-Co chains, the Co-0 bond 
is about 0.18 A longer than the Mn-0 bond. This clearly 
indicates that the Mn-0 bond is stronger in comparison 
with Co-0 bond, which is reflected in the DOS where we 
see that more Mn-d states are pushed down below the 
Fermi level. 

We have computed a spontaneous electronic polariza- 
tion of 5.88 [iCjcw? in the FM state by Berry phase 
method. A 4 x 4 array of strings containing 32 k-points 
has been used to obtain the electronic polarization, which 
is non-zero only along the short-axis b in agreement with 
the C2 symmetry of the crystal. Our estimated value is 
close to what has been measured for Bi2NiMn06 thin 
films. Unfortunately, we are unable to find an experi- 
mental value for BCMO. 

Now we discuss the effect of magnetic structure on the 
BEG matrices shown in Table 1. It is observed that they 
are quite anisotropic due to the highly distorted structure 
of the compound and there are significant differences in 
BEC's for FM and FiM states for Co and Mn. Both di- 
agonal and off-diagonal components of BEG matrices of 
Co and Mn are different in FM and FiM states, with ma- 
jor changes observed for [Z*')j,\ and {Z*)32 components 
(shown in box) of Co and Mn BEC's respectively. These 
changes may be attributed to the changes in electronic 
DOS's in two magnetic states, where, e.g., the energy 
separation between the occupied Co and unoccupied Mn 
states in the minority spin channel (Fig. 2) is de- 
creased significantly in FiM state compared to the FM 
state. According to Ghosez et al^, the expression for 
BEG has the energy separation in the denominator and 
hence a reduction of this value should increase BEG sig- 
nificantly. A relatively small change occurs in the di- 
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agonal components of Mn as the Mn-0 bonding is al- 
most uniform over a wide range of energy whereas a sig- 
nificant Co-0 hybridization about 0.5 below the Fermi 
energy increases the diagonal component of Co EEC's. 
The EEC's of O are almost identical in both structures. 
Among all the species, Ei has the largest EEC, demon- 
strating a strong Ei-0 covalancy. From all these observa- 
tions we can conclude a significant magnetism dependent 
electrostructural coupling in the system very similar to 
what has been observed by Das et. aP^for La2NiMn06. 
However unlike the study in Ref. 1161 we did not observe 
an anomalously large off-diagonal EEC for either of the 
magnetic specie. 
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FIG. 3: Phonon DOS of IR-active modes for FM and FiM 
configurations. 

Finally, we analyze the zone center phonon modes for 
FM and FiM structures. The eight lowest infra-red (IR) 



active modes are shown in Fig. 3. One can observe that 
the typical softening of IR-active modes could be as large 
as 10 cm~\ The frequencies at 62, 66, 79, 82, 87, 133,133 
and 169 cm^^ in the FM phase are shifted to 52, 54, 61, 
67, 83, 124, 129 and 156 cm^^ in the FiM phase. The 
least shifted mode is at 133 cm^^ where the softening is 
associated with the further lifting of degeneracy of this 
mode. Subsequently there is an increase in the static 
dielectric constant in FiM phase. To understand the spin- 
phonon coupling better, we analyzed the softest mode (at 
52 cm"'^) in the FiM state more carefully. Ey adding the 
atomic displacement due to this mode to the equilibrium 
atomic positions we find that the atomic displacements 
are such that the increase in Mn-O-Co angle with respect 
to FM state is very small. We find an average increase 
of about 2° of Mn-O-Co angle. This fact may contribute 
to the observation that the spin-phonon coupling though 
present, is small. Moreover, the contribution of Mn to the 
change in mode effective charges of the IR-active modes 
of two magnetic states as well as to dielectric constant is 
smaller compared to Co. 

In summary, our GGA-I-U calculations reveal that Co 
and Mn are in high and intermediate spin states re- 
spectively in the double perovskitc Ei2CoMn06 with a 
calculated band gap of 0.3 eV. The computed dynam- 
ical charge matrices are highly anisotropic due to the 
low symmetry structure. Magnetic structure dependent 
phonon frequencies indicate a spin-phonon coupling in 
this compound. Our theoretically estimated spontaneous 
polarization of 5.88 /cm? is close to the experimental 
value observed for the thin film of a similar compound, 
EiaNiMnOe. 
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